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Abstract 

The coupled Einstein- Yang- Mills equations on a time dependent axially sym- 
metric spacetime are investigated, without imposing a priori any conditions 
on the gauge field. There is numerical evidence for the existence of a regu- 
lar solution with the desired asymptotic features. Just as in the supermassive 
abelian counterpart model, the formation of a singularity at finite distance of 
the core of the string depends critically on a parameter of the model, i.e., the 
constant value of one of the magnetic components of the YM potentials. The 
multiple-scale method could supply decisive answers concerning the stability 
of the solution. 
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I. INTRODUCTION 



After the astonishing numerical solution of Bartnik and McKinnon of the static spher- 
ically symmetric Einstein- Yang-Mills equations with a SU(2) gauge group, a fundamental 
question arose: do there exist essentially non-abelian static globally regular solutions or 
blackholes, which carry non-zero electric and/or magnetic charges In first instance, 

one could prove in the spherically symmetric case |^ the existence of static solutions which 
are asymptotically flat and with finite mass. The members of this family of solutions could 
be characterized by the number of zeros of one of the gauge potentials. In fact the repulsive 
Yang-Mills force can balance the gravitational attractive force and prevent the formation of 
singularities in spacetime. Later p the class of solutions were extended to blackhole solu- 
tions with a horizon at some and oscillating solutions which are not asymptotically fiat. 
These non-abelian blackholes could represent counter examples to the 'no-hair' conjecture. 
Further, it was claimed that both regular and blackhole solution in the spherically sym- 
metric case will have gravitational-like and sphaleron-like instabilities. When it was found 
that the BM solution is unstable, a lot of researchers (see for example tackled the 

general stability problem of the coupled EYM system and investigated the possible critical 
behavior of the solutions. It was realized that they had much in common with the elec- 
troweak sphalerons. In the spherical symmetric case the critical behavior was formulated 
by Choptuik [§,0]. It occurred at the boundary in phase space between initial data which 
eventually form a blackhole and data which do not. Recently it was found, in a two dimen- 
sional Weinberg-Salam model with an axially symmetric ansatz, that there is evidence for 



the existence of an electrically charged sphaleron state [Til]. H^re we consider the EYM sys- 
tem on an axially symmetric time dependent spacetime and will compare our investigations 



with the spherical symmetric case and the well known abelian cosmic string solution |12-15 



particularly the supermassive case. In these U(l)-gauge cosmic string models it was found 
that as the energy scale of symmetry breaking increases, the geometry around the string 
changes from conical to an analog of a Kasner spacetime. But supermassive cosmic strings 
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may also arise at GUT scales if the coupling between scalar and gauge fields is weak. Fur- 
ther, these low-energy supermassive strings are closely related to global strings: they both 
show singular behavior at finite distance from the core of the string. We conjecture that 
the string-like solution in the Eintein- Yang-Mills model shows some similar behavior. The 
plan of this paper is as follows. In section II we derive the time dependent equations for 
the Einstein- Yang-Mills system on an axially symmetric spacetime using the algebraic ma- 
nipulation program MAPLE. In section III we consider the static situation in order to gain 
insight into the asymptotic and singular behavior. In section IV we solve the equations 
numerically and in section V we present a conclusion and prospect concerning the stability. 
We conjecture that conventional (linear) stability analysis is inadequate to be applied to the 
situation where a singularity is formed. The reason is that the oscillatory behavior of some 
of the gauge field components are of high-frequency. We suggest to apply the multiple-scale 



method |16-18|, a method suitable for handling high-frequency perturbations. In a future 



work, we will extend the investigation to the spinning case, as initiated before |p!9| , p0 | 



II. THE FIELD EQUATIONS 



Consider the Lagrangian of the SU(2) EYM system 

TZ 1 

167rG ~ 4' 



S = I (fx^-detg - -T^T^""" , (2.1) 



with the YM field strength 

T% = d,Al - d.Al + ge'^^'AlAl, (2.2) 

g the gauge coupling constant, G Newton's constant, A" the gauge potential and 7?. the 
curvature scalar. The field equations then become 

G^, = -SttGT^,, (2.3) 

^^JTf^ua ^ ^2.4) 
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with T the energy momentum tensor 



1 



= ^U^l"" - -^9,.Kp^''^^- (2.5) 
We are interested in solutions on the cyhndrical symmetric spacetime 

ds" = -e2(^-^)(cit2 _ dr'') + e-^^'W^d^^ + e^"" dz\ (2.6) 

where t/, K and W are functions of t and r. This cyhndrical symmetric line element can 
formally be obtained from the stationary axisymmetric line element used by Kleihaus and 
Kunz m 

ds^ = -fdt' + j{dr' + dz') + y V, (2.7) 



with /, m and / functions of r and z, by the complex substitution t iz and z it ^ 



Following Bais and Sasaki we then can specify the gauge field potentials = A^t^, 



where the ta are the anti-hermitian generators of the gauge group, by 

A^dx^ = Aof^dt + Aif^dr + [7727V + r]if^]dz + [x2Tr + Xi^z]d'^, (2.8) 

with Ai^rji and Xi functions of t and r and where f^,fr and are the axial generators 
of the SU(2) normalized such that [ri^fj] = eijufk- One can reduce the number of gauge 
field potentials. On the Euclidean spacetime, the self-duality condition eliminates two of 



the six functions |24]. One also can use the additional U(l) gauge freedom on the A^, 
A — i> h^^Ah + h^^dh, with h = exp['0(r, t)f^] 0. Here we use the condition that on the 
spacetime (2.6), the energy-momentum component T^ip must vanish. The most simple way 
to fulfil this condition is rji = X2 and 7/2 = Xi- There are some other possibilities, comparable 
with those found by Manton P||. However, these are quite complicated as already noticed 



by Manton. So in our case, T^^p = becomes 

%^ = gr]i{Al - Al){l + 2gr]2) + 2(a,r/i9,772 - dtr]idtr]2) + Aodt7]2 + AidrV2 = 0. (2.9) 

We then have two possibilities: i. 772 =constant= 7/0 and Aq = ±Ai, or, ii. V2 = We 
consider here case i. The condition of case ii means that the YM mass scale Mym = gVo = 
— ^. In section 3 we shall see that it corresponds in the static case with a singular solution. 
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The differential equations for the several field variables become (from now on we set 
Vi = v) 

d^T] - dlf] = ^dr'n^drW - 2WdrU) + l^dtvi2WdtU - dtW) - gvoiOrAo + dtAo)+ 

A '2K 

^"^^ '{2WdtU + 2WdrU - dtW - drW) + ^^[m + givo - v% (2.10) 



w 



^2K 



O T T O T T 7" 

+{dtAo + 9,Ao)(^- + 2drU - 2dtU + 2dtK - 2drK) 



-2 



W W 

"^[gW + (1 + griof] - 2e-^-^^ A,g-{jf + ryo^), 



(2.11) 



dlW - 



SttG 



W 



^'^-'^b(%' - V') + ^0]^ + e-™W^^(9,Ao + d^AoY , (2.12) 



d'^K - d^K = {dtUf - {drUf + 



47rG 



-2C7+2Kr„/ 2 _ 2 



+2^obr;o^'e-2^ + {gi^o + l)e'^](a,77 + dtvi) 



(2.13) 



and 



dlU - dJU = dtUdt InW- dMdr InW - 



AttG 



-2U+2K:„(J2 _ -1 



-2Ao[gvoW^e-^'' + (grjo + l)e^'']{drV + dtv) • 
Prom two combinations of the Yang-Mills equations (2.4), i.e., 

[yM],=^,„=3 ± (cosv?.[yM],=,,„=i + sin(^.[yM],=,,„=2), 

we obtain an expression for ^40, 

_ e^^-^^W'ivl + gyoivl - ry^)] + e^^[ai(l + r^^gf - r,\l + grjo)] 
° ~ 2r)W[dtW + drW - 2W{drU + dtU)] 

and a first order equation for Aq 



(2.14) 



(2.15) 
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dtAo + drAo 



Tj Irj 



■[gr] {l + gr]o) - r]o{l + gr]o) ]. 



(2.16) 



2gr]W^ 

In the abelian situation, the equations could be simphfied, due to the condition Tu = —T^z-, 
and one obtains [TH] K = 2U. This is the famihar self-gravitating Nielsen-Olesen vortex 



model studied first by Garfinkle [|1J] in the static case. In our case we obtain the condition 



2-2K 



(2.17) 



which yields 

dlK - dlK = dlu - d^u + ^{drUdrW - dtUdtW) - {d^uf + {dtUf - ^{olw - d^w). 

(2.18) 



This relation can be used as constraint in the numerical code. 



III. THE STATIC CASE 



In order to obtain boundary conditions for the numerical integration of the differential 
equations, we first consider the static situation. The field equations reduce to 
-SttG 



2 J2\i2 



(3.1) 



dr[WdrU] 



W 



(V[/2e-2C/ _ e2f^)(9,,?7)2 + W^e'^^'-^^'idrAo) 



r2-2U 



2K-2U\ 



.2 „2m2 



(3.2) 



dr[WdrK] 



8nG 

~w 



{g\ri + r/^)A^l^^e-2^ + + (1 + g^,)A,Y + g^e'^^'K , (3.3) 



2 I „2 2V 2 a2 



--drAo + 



[1 + 2gr]o)Ac 
2W 



{drW - 2WdrU) 



2W^ 



(3.4) 



and 
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dlAo = drAo{2drK - ^ - 2drU) + ^e^^[(l + gvoY + gW] 



Further, (2.15) now reads 

e'^-'^W'[v', + grioivi - y')] + e^^[f (1 + gyof - v\l + QVo)] 
° ~ 2r]W{drW - 2WdrU) 

From the condition that in the static situation now Ttr = 0, we obtain 

which is consistent with one of the YM equations. For the special case gr]^ = — | 
from (3.7) 



(3.5) 



(3.6) 



drT] 



^)2 + l 2 



1 , 



(3.7) 
we obtain 

(3.8) 



which can be integrated together with (3.6) for some given spacetime. However, we observe 
from (3.6) that Aq becomes singular for W — > e^^ . We then obtain from the YM equations 
just the condition grjQ = —\- 



The asymptotic features of the system can be analyzed by following Garfinkle [|T^]. We write 
the stress tensor on an orthonormal basis as 



where t^' = e-^+^{^Y,f^ = e-^+^(|:)^,5 
vation of stress energy 



(3.9) 



l'^. Using the conser- 



v.r^" = 0, 



(3.10) 



one obtains 



d_ 

dr 



WVr ={a + Vr + V,- V^)Qi - (a + P,)03 + ^^©2, 



(3.11) 



where 
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Gi = WdrU, 82 = drW, 83 = WdrK. (3.12) 

Let us assume that 



and that 6j approach constant values as r — > 00 ||T2|, which are fairly weak assumptions. 
Then one can write with the help of (3.1)-(3.3) 

One then obtains the asymptotic condition 

lim (6263 - e?) = 0. (3.15) 

r— >oo 

So one has the two possibilities 

i : WdrU\^ = 0, WdrWdrK\^ = 

ii: {WdrUlooY = WdrWdrK)\^. (3.16) 

If we denote by W^o , and i^oo the values of the metric fields far from the string, one 
obtains in the first case: Uoo = fli, Woo = a2r + a3 and Koo = 0,4, where the are constants. 



The metric approaches in this conical spacetime lO 



ds^ = e2('^4-ai)^_^^2 ^ ^^2] ^ g2ai^^2 ^ ^~2a^ ^^^^ ^ tt^fd^^ (3.17) 



In the second case one obtains in the special case K = 2U the Kasner-like spacetime 

p-263 

ds^ = -{h^r + h^fe'^'^ide - dr^ - dz^] + -———d^\ (3.18) 

(6ir + 62) 

with hi constants. Now we can investigate the asymptotic behavior of 7] and Aq on the conical 
spacetime (3.17) and compare the asymptotic features with the abelian counterpart model. 
Substituting the conical spacetime (3.17) into (3.6) and (3.7) one obtains solutions for rj and 
Aq. The solution is not sensitive for ai, 03 and 04. For positive 02 the solution will possess 
singular behavior at finite distance of the core, as expected from the abelian investigations 
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of Laguna and Garfinkle [|T^]. In figure 1-3 we plotted rj and Aq for CLl — CI2 — O'A — 

1 and 

03 = 2 and for different values of tjq. It is observed that the singularity is pushed to infinity 
for smaller values of rjQ and 02 > 1. This kind of singular behavior is also encountered in 



the abelian Higgs model ||T3|,|T^ when the scale of the symmetry breaking is far beyond the 
GUT scale, i.e., SttGi^^ >> 10^^, where u is the vacuum expectation value of the Higgs field. 
In fact, the conical picture of the string then fades away. It was found that the metric of 
the transition between conical and Kasner-like is cylindrical, i.e., i?^ x S*^. Further it was 
found that supermassive strings can be formed at GUT scales if a = ^^^^^ , where e is the 
gauge coupling constant in the abelian model and A the Higgs self coupling constant, is very 
small. The singularity must then occur at a distance from the core of the string which is 
many orders of magnitude greater than the present Hubble scale. In our model we find a 
comparable dependency on tjq of the behavior of the singularity. In general, however, one 
has to solve simultaneously the coupled system. Moreover, one should make a distinction 
between the interior and exterior field equations with proper matching conditions, in order 



to get insight in the behavior close to the z-axis ||20|| . 



IV. NUMERICAL SOLUTION 

The set of equation (2.10)-(2.14) can be solved numerically. We used the ISML software 
package for numerically solving coupled systems of nonlinear partial differential equations. 
The package implements finite element collocation methods based on piecewise polynomials 
for the spatial discretization techniques. The time integration process is then accomplished 
for a set of ordinary differential equations using banded Jacobians. For the order of the 
piecewise polynomial space we took 5 and for the number of subintervals into which the 
spatial domain is to be divided we took 11. The relative error bound was 10^^^. In order 
to obtain from (2.15) a regular and asymptotically correct initial value for Ao, we choose as 
initial values r]{r, 0) = e"^ + 1, U{r, 0) = K{r, 0) = and W{r, 0) = r + 1. Fi gures 4, 5 and 
6 show a typical regular solution of Aq,!] and the metric component e^^~^^ for some values 
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of r]o and AirG and suitable boundary conditions. The behavior of 77 remains regular every- 
where. In figure 7 we plotted Aq for a smaller value of tiq. We observe that Aq approaches 
asymptotically a constant value. For runs where Aq starts to oscillate, e^^'"^^ decreases 
strongly, signalling the formation of a singularity. However, a thoroughly investigation of 
the dependency on rjo and G will be necessary. This is currently under study. 



V. CONCLUSIONS AND OUTLOOK 

We investigated the existence of a possible solution of the coupled EYM system on a time- 
dependent axially symmetric spacetime. At least on the conical spacetime we find evidence 
for regular behavior of the electric and magnetic components of the YM field for suitable 
values of the parameters 'qQ, g and AttG. The formation of a singularity at finite distance of 
the core of the string depends critically on r]o, the constant value of one of the gauge field 
potentials, and AnG. Just as in the static spherically symmetric EYMH case, where the 
ratio = is the vacuum expectation value of the Higgs field) plays a crucial role 

in the behavior of the solution, we find a crucial dependency on the ratio ^^J^' = J^°q - In 



the abelian supermassive model a similar behavior is encountered where the singularity 
arises at finite distance from the core of the string not only for large symmetry breaking 
scales as found by |13|, but also for GUT scale. There is for fixed symmetry breaking scale of 



order of the GUT scale a critical value for the coupling of the scalar to the gauge field (i.e., 
^^^'^^ with e the coupling constant and A the Higgs self coupling) for which the singularity 
occurs at finite distance of the core. These low energy solutions are more realistic because 
the larger the energy scale the larger the angle deficit and it no longer makes sense to talk 
about string type solutions. In our model, there will be a critical value J^^^ for which the 
singularity is pushed to infinity. 

In order to analyze the stability of the solutions, one usually linearizes the field equa- 



tions B, or one expands the field variables in a physically unclear small parameter |25]. One 



better can apply the so called multiple-scale (or two-timing) method, developed decades ago 
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by Taub |T6[ and Choquet-Bruhat [H^. This method is is particularly useful for constructing 



uniformly valid approximations to solutions of perturbation problems. The idea is to expand 
the several field variables in power series of the ratio of the characteristic wavelength of the 



perturbations and the characteristic dimension of the background. One writes 18 



Al = Al + + + (5.1) 



where ^ = ull{x'^) and 11 a phase function. The parameter u measures the ratio of the fast 
scale to the slow one. The rapid variation only occur in the direction of the vector 1^ ~ 



dx" 



For a function '^{x"; C,) one has 

dx"" 



d^'^ + uL^, (5.2) 



where d^"^ = -§^\^fixed and \i/ = ^Ix" fixed- Substituting the expansions of the field vari- 
ables into the equations and collecting terms of equal orders of uj, one obtains propagation 
equations for B'^ and h^^, and 'back-reaction' equations for h^i_, and It is clear from the 
propagation equation that there will be a coupling between the high-frequency gravitational 
field and the high-frequency behavior of Aq when the singularity will be approached. On his 
turn, Aq will create a high-frequency perturbation in rj. 
In a subsequent paper we will present this investigation. 
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FIGURES 



FIG. 1. Plot of the gauge components Aq and rj for rjo = 0.5. The singularity appears when 77 
approaches zero 
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FIG. 2. As figure 1, with 770 = 0.2. The singularity is encountered at larger r value 



15 



FIG. 3. As figure 1, with ijq = 0.01. The singularity is pushed to infinity 
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FIG. 4. Plot of a long time-run of the gauge component AO for g=-l, 'iirG = 0.2 and rjo = 0.2 



FIG. 5. Plot of ij for the situation of figure 4 
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FIG. 6. Plot of the metric component e for the situation of figure 4 



FIG. 7. As figure 4, but now for rjo = 0.01 
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